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Background: The respiratory tract is a major target of exposure to air pollutants, and respiratory 
diseases are associated with both short- and long-term exposures. We hypothesized that improved 
air quality in North Carolina was associated with reduced rates of death from respiratory diseases 
in local populations. 

Materials and methods: We analyzed the trends of emphysema, asthma, and pneumonia 
mortality and changes of the levels of ozone, sulfiir dioxide (SO^), nitrogen dioxide (NO^), carbon 
monoxide (CO), and particulate matters (PM^^ and PM^^) using monthly data measurements 
from air-monitoring stations in North Carolina in 1 993-20 10. The log-linear model was used to 
evaluate associations between air-pollutant levels and age-adjusted death rates (per 100,000 of 
population) calculated for 5-year age-groups and for standard 2000 North Carolina population. 
The studied associations were adjusted by age group-specific smoking prevalence and seasonal 
fluctuations of disease-specific respiratory deaths. 

Results: Decline in emphysema deaths was associated with decreasing levels of SO^ and CO in 
the air, decline in asthma deaths-with lower SO^, CO, and PMjjj levels, and decline in pneumonia 
deaths-with lower levels of SO^. Sensitivity analyses were performed to study potential effects 
of the change from International Classification of Diseases (ICD)-9 to ICD-10 codes, the effects 
of air pollutants on mortality during summer and winter, the impact of approach when only the 
underlying causes of deaths were used, and when mortality and air-quality data were analyzed 
on the county level. In each case, the results of sensitivity analyses demonstrated stability. The 
importance of analysis of pneumonia as an underlying cause of death was also highlighted. 
Conclusion: Significant associations were observed between decreasing death rates of emphysema, 
asthma, and pneumonia and decreases in levels of ambient air pollutants in North Carolina. 
Keywords: chronic obstructive pulmonary disease, sulfur dioxide, carbon monoxide, nitrogen 
dioxide, particulate matter 

Introduction 

Air pollution has a deleterious impact on human health,'"' with global outdoor air 
pollutants estimated to account for approximately 1 .4% of total mortality and 2% of 
all cardiopulmonary mortality.^ Both ambient particles'*'* ' and such gases as nitrogen 
dioxide (NO^), ozone (Oj), and carbon monoxide (CO) have been shown to increase 
total, cardiovascular, and respiratory (predominantly due to lung cancer and chronic 
obstructive pulmonary disease [COPD]) mortality and morbidity.^ '" " While the impact 
on any individual's risk of death has been thought to be relatively modest per se, the 
overall impact of air pollution on the health of an exposed population makes it a major 
public health concern. 
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While more studies on short-term impacts of changes 
of air quality are available (such as the legislated traffic 
holidays during the 1996 Atlanta Olympic Games'^ and the 
2008 Beijing Olympic Games'''), less is known about the 
long-term effects of changing air quality on the health of 
exposed populations. For example, a ban on heating-coal 
sales in Dublin was thought to be associated with both 
reduced pollution from airborne particulate matters (PMs) 
and 5.7% reduction in all-cause, 15.5% reduction in respi- 
ratory, and 10.3% reduction in cardiovascular mortality.'^ 
However, these results were considered inconclusive, due to 
the complexity and expense of evaluating the health effects 
of air pollution on populations."" '^ Since the 1990s, a variety 
of acts, standards, and requirements in the US have been 
adopted to improve air quality. For example, increasingly 
stringent national gasoline and automotive engine require- 
ments have been applied, resulting in a decrease of CO, NO^, 
PM, and volatile organic compounds in the air. At the state 
level. North Carolina in 1 992 entered into the Southern Appa- 
lachian Mountains Initiative, leading to the development of 
the Clean Smokestacks Act'* to mandate reduced emissions 
from coal-fired power plants." 

While few studies have analyzed the associations of 
both air quality and health over a long period, and they were 
typically limited to analysis of a specific air pollutant or a 
couple of pollutants, we were able to study longitudinally a 
number of air contaminants, including both PMs and noxious 
gases. In addition, we analyzed both air quality and health 
outcomes over almost two decades (1993-2010). Because 
respiratory morbidity and mortality are affected by changes 
in air quality,^"'^^ we evaluated the associations between the 
changes of the levels of PMj^ and PM^ 5, ozone, CO, NO^, 
and SOj in the air and death rates of emphysema, asthma, 
and pneumonia. 

Materials and methods 

Data 

We analyzed mortality rates for emphysema {International 
Classification of Diseases [ICDJ-9 code 492, ICD-10 code 
J43), asthma (ICD-9 code 493, ICD-10 codes J45, J46), and 
pneumonia (ICD-9 codes 480.0, 480. 1 , 480.2, 480.9, 485, 486, 
487.0, 487.1, ICD-10 codes Jl 1.00, Jl 1.1, J12.0, J12.1, J12.2, 
J12.9, J18.0, J18.9) inNorth Carolina from 1983 to 2010 using 
the data fi-om the Vital Statistics National Center for Health Sta- 
tistics Multiple Cause of Death dataset. We started the mortality 
analysis with the data from 1983, but could only analyze air 
quality when monitoring data were available, ie, 1993-2010. 
The mortality data enabled an analysis of a longer period of 



death-rate dynamics, thus allowing to observe the dynamics of 
disease-specific mortality before the measured reduction in par- 
ticulate and gaseous emissions in North Carolina. Age-adjusted 
death rates (per 100,000 of population) were calculated using 
5-year age-groups and standard 2000 North Carolina popula- 
tion. The data on population were provided by the Surveillance 
Epidemiology and End Results Registry (SEER) at http://www. 
seer.cancer.gov/popdata/download.htmL 

Data on concentrations of PM^^ (|j.g/m^), PMj|, (|j,g/m^), 
ozone (ppb), CO (ppb), NOj (ppb), and SO^ (ppb) in the air 
in 1993-2010 were obtained from the US Environmental 
Protection Agency (EPA) ( http://www.epa.gov/ttn/airs/air- 
saqs/detaildata/downloadaqsdata.htm ). We used the averaged 
month-specific concentrations of air pollutants for North 
Carolina to further analyze them for associations with the 
dynamics of cause-specific monthly mortality in the state. 
A two-stage averaging procedure was used to avoid heteroge- 
neity in the numbers of measurements made in certain days 
of the month: first, we calculated the day-specific means, and 
then these values were averaged, resulting in month-specific 
means. Negative values were excluded, and measurements 
with various units were converted to |J,g/m' for PM^^ and 
PMji,, and to ppb for ozone, CO, NO^, and SO^. Since the data 
on air pollutants represented different methods of registration 
during different durations of sample collection (ie, the length 
of time used to acquire a sample measurement), an auxiliary 
analysis was performed to check whether the specific method 
could be considered as an outlier and therefore excluded 
from the analyses. 

Also, data on the prevalence of tobacco use for 1 995-20 1 0 
were obtained from the Centers for Disease Control and 
Prevention Behavioral Risk Factor Surveillance System 
survey for age-groups 18-24, 25-34, 35^4, 45-54, 55-64, 
and 65+ years ( http://www.cdc.gov/brfss ). 

Ethics statement 

The data used in this study have no individual identifiable 
information. No specific procedures were required for de- 
identification of the records. All data analyses were designed 
and performed in accordance with the ethical standards 
of the committee on human experimentation and with the 
Helsinki Declaration (1975, revised in 1983), and were 
approved by the Duke University Health System Institutional 
Review Board. 

Methods 

Trends of cause-specific death rates and of levels of air 
contaminants were analyzed for correlations. Adjustment by 
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smoking prevalence and seasonal fluctuations in respiratory 
deaths (for monthly death rates of emphysema, asthma, and 
pneumonia) were included in a log-linear model that was 
used to evaluate the associations between the level of each 
studied air pollutant and the death rates, as follows: 

11 

log(r) = M + ^c + ^2^ + ^//y„+f, (1) 

m=l 

where u was the intercept, represented the effect of each 
studied air pollutant depending on its concentration (denoted 
by c) measured in its units (as described in the Data section), 
represented the effect of smoking prevalence (denoted 
by s), represented the effects of 1 1 months (January to 
November for each year) in respect of December (/^ is the 
month indicator), and e stood for random residuals. Note 
that if the air-pollutant concentration changes by one unit of 
its measured level in the air, the rate r changes by the factor 
of exp(/?j). For multiple comparisons, the Bonferroni cor- 
rection was applied. 

Sensitivity analysis 

The potential effect of ICD code changes (from ICD-9 to 
ICD-10), the seasonal fluctuation of air pollutants and mor- 
tality during summer and winter, and the analysis validity 
when only the underlying causes of deaths contributed to the 
cause-specific death rates were tested. In addition, sensitivity 
analysis was performed for county-level data on respiratory 
mortality and air-pollutant levels. Only counties for which 
the data on air quality were directly measured by monitoring 
stations were included in the analysis: 37 counties for ozone 
measurements, 1 1 counties for NO^, 22 counties for SO^, 16 
counties for CO, and 37 counties for PM^ ^ and PMju measure- 
ments. As in the main analysis, dynamics of smoking preva- 
lence (on state level) and seasonal fluctuations in respiratory 
mortality were used for adjustments of the results. 

Results 

We analyzed up to 1 80 month-specific measurements of each 
of the studied air pollutants recorded at multiple monitoring 
sites in North Carolina (see Table 1 for detailed air pollutant- 
specific information). We found air quality in North Carolina 
gradually improving over time, primarily due to decreasing 
PMjj, NOj, and CO levels. These decreases became more 
pronounced from 2002 (see Figure 1 ; note that individual 
pollutants were placed onto a single graph by utilizing the 
arbitrary units to enable a collective visualization of the 
trends). The following seasonal fluctuations of pollutants 
levels were observed (Figure 2): levels of ozone, PM^^, and 



Table I Measurements of air pollutants used in the study, 
1993-2010 



Air pollutant 


Number of 


Number of month- 




monitored sites 


specific measurements 


Ozone 


69 


148 


Nitrogen dioxide 


IS 


180 


Sulfur dioxide 


3S 


180 


Carbon monoxide 


41 


180 




68 


180 




60 


132 



Abbreviation: PM, particulate matter. 



PMji, were higher in summer, while levels of SO^, NO^, and 
CO were higher in winter. 

Since 1 983, the death rates of three studied diseases have 
been decreasing (Figure 3), with declines in emphysema 
death rates more dramatic since 1998, for asthma since 1995, 
and for pneumonia since 1990. From 1993 to 2010, 101,374 
deaths in North Carolina were caused by pneumonia, 13,187 
by emphysema, and 5,509 by asthma. The detailed description 
of the studied population is presented in Table 2. Among those 
who died from emphysema and from pneumonia, 80.7% and 
85.9%, respectively, were older than 65 years. For asthma, 
ages at death were younger: 9.7% were younger than 40 years, 
and 3 1 .3%) were aged 40-64 years old. However, the declin- 
ing trends of pollutant concentrations and death rates during 
1993-2010 do not essentially confirm causality. 

The association between the changes of air-pollutant 
levels and dynamics of disease-specific death rates after 
being adjusted for smoking prevalence (for respective year 
and age-group), and by monthly fluctuations in respiratory 
disease-specific death rates are shown in Table 3, for each 
air pollutant. The disease-specific death rate (number of 
deaths per 100,000 population) decreased by a factor cal- 
culated based on the value of estimate presented in Table 3 
(ie, per decrease of concentration of each pollutant by one 
unit of measurement: per 1.0 ppb for ozone, SO^, NO^, CO, 
and per 1.0 |J,g/m^ for PM^, and PM^^). For example, the 
estimate for emphysema in Table 3 means that if the SOj 
level decreases by 1 ppb, the emphysema death rate (per 
100,000 population) can be predicted to decrease by a factor 
of exp(0.0547) =1.056. Similar interpretation can be devel- 
oped for smoking estimates, keeping in mind that smoking 
is represented by its prevalence in population measured in 
percentages, and thus the respective exponential factor cor- 
responds to a change in smoking prevalence by \%. 

Among gaseous pollutants, the estimates for associa- 
tions between reduction of air-pollutant levels and reduc- 
tion of death rates were significant for SO^ and emphysema 
(0.0547+0.0106, P<0.0001), asthma (0.0598+0.0173, 
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Figure I Levels of six air pollutants in North Carolina, 1993-201 I. Individual pollutants were placed onto a single graph by utilizing arbitrary units to enable a collective 
visualization of the trends. 
Abbreviation: PM, particulate matter. 



/'<0.001), and pneumonia (0.0309+0.0093, P<0.001), and 
for CO and emphysema (0.0004+0.0001, P<0.0001) and 
asthma (0.0006+0.0001, P<0.001). For PM, reduced PM^^ 
levels were associated with reduction of emphysema mortal- 
ity (0.0155+0.0066, /'<0.05) and reduced PMj„ levels, with 
reduction of asthma mortality (0.0204+0.0058, /'<0.001). 
As expected, smoking significantly affected the mortality of 
each disease. 

Sensitivity analysis 

The sensitivity analysis demonstrated good stability of obtained 
results (see Table S 1 for detailed information). In the sensitivity 
analysis, the association between pneumonia mortality and CO 
levels became significant (/'=0.0655 in main versus P<0.0001 
in sensitivity analysis) when pneumonia was analyzed as an 



underlying cause of death. Recent studies have demonstrated 
that separation of comorbid conditions to underlying and 
secondary causes can be unreliable;^'"^^ however, for certain 
diseases with a predominantly acute course (eg, pneumonia), 
that may not be the case, and additional information can also 
be obtained from analysis of underlying causes of death. In 
addition, sensitivity analysis showed that during summer 
decreased mortality from emphysema was associated with 
lower levels of PM^g (P=0.2554 in main versus P=0.017 in 
sensitivity analysis), and statistical significance was observed 
for associations between pneumonia mortality and CO levels 
when ICD code changes were taken into account (P=0.0655 
in main versus P=0.018 in sensitivity analysis). 

A county-level analysis also demonstrated the stability of 
most observations in the main analysis. Among associations 
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Figure 2 Seasonal fluctuations of air-pollutant levels: summer (red, 3 months) and winter {blue, 3 months), 1993-201 I. 
Abbreviation: PM. particulate matter. 
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Figure 3 Trends in death rates for emphysema, asthma, and pneumonia in North Carolina, 1983-2010. Mortality rates were age-adjusted to the 2000 North Carolina population. 
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Table 2 Demographic characteristics of North Carolina 
population with cause-specific respiratory mortality, 1993-2010 



Demographic 
characteristic 


Cause of death 




Emphysema 


Asthma 


Pneumonia 


Number of deaths 


13,187 


5,509 


10,1374 


Sex, n 








Males 


7,951 (60.3%) 


1,806 (32.8%) 


48,517 (47.9%) 


Females 


5,236 (39.7%) 


3,703 (67.2%) 


52,857 (52.1%) 


Race, n 








Caucasians 


1 1,866 (90.0%) 


3,567 (64.8%) 


82,759 (81.6%) 


African-Americans 


1,237 (9.4%) 


1,853 (33.6%) 


17,665 (17.4%) 


Other 


84 (0.6%) 


89 (1.6%) 


950 (1.0%) 


Age, n 








< 1 5 years old 


7 (0.1%) 


103 (1.9%) 


595 (0.6%) 


1 S-39 years old 


40 (0.3%) 


429 (7.8%) 


1,633 (1.6%) 


40-64 years old 


2,504 (19.0%) 


1,723 (31.3%) 


12,054 (1 1.9%) 


65+ years old 


10,636 (80.7%) 


3,254 (59.1%) 


87,091 (85.9%) 



that were significant under Bonferroni correction in the main 
analysis, associations between dynamics of SO^ and mor- 
tality from emphysema (0. 1399, P<0.00\) and pneumonia 
(0.0698, /'<0.001), and associations between changes of 
CO levels and asthma mortality (0.0004, P<0.05) were 
also significant in the sensitivity analysis. The association 
between CO and pneumonia mortality was also significant 
when analysis was performed on a county level (0.0002, 
P<0.00\). Recall that this association was significant in 
the analysis using state-level data in two cases: when being 
corrected for changes of ICD codes and when only underly- 
ing causes of deaths were considered as contributing to the 
cause-specific death (see detailed results in Table SI). The 
effects of dynamics of SO^ and PM^^ on asthma mortality 
became nonsignificant {P>0.05), likely due to the small 
number of county-specific asthma deaths and due to the 



large fraction of zeroth death rates that were not successfully 
described by Equation 1 . 

Discussion 

We foimd significant correlations between reduction of air 
pollutants and dynamics of deaths due to respiratory diseases 
during the period we studied. We need to contextualize 
these findings, particularly in regard to the multifactorial 
contributors to respiratory mortality. In general, COPD has 
been shown to correlate highly with air pollution linked 
to global urbanization,^' eg, higher prevalence of chronic 
bronchitis (odds ratio [OR] 2.26, confidence interval [CI] 
1 .54-3.3 1), asthma (OR 1 .57, CI 1 .25-1.98), and emphysema 
(OR 2.98, CI 1 .95^.54) were observed in the meta-analyses 
of individuals exposed to urban air.^' Little is known about 
whether chronic, low-dose exposure to ambient air pollutants 
can exacerbate COPD progression.^*'^' Several recent stud- 
ies related respiratory symptoms to long-term rather than 
short-term effects of ambient particles,'" with the long-term 
exposure to PM^^ increasing the risk of COPD.'' 

Changing air quality in North Carolina could be a good 
example of analysis of the trends of both improved air quality 
and respiratory mortality over almost two decades of obser- 
vations. Improved air quality in North Carolina since the 
mid-1990s is related to a series of federal and state acts and 
regulations (see Table 4), including the national heavy-duty 
truck engine standards, reduction of NO^ emissions, the Clean 
Smokestacks Act, and new engine standards. Regulations of 
emissions of NO^, PMj|,, and CO appeared to be very effec- 
tive in improving air quality in the state. Observed seasonal 
fluctuations of air-pollutants levels could be due to season- 
dependent local dispersive conditions, breeze dynamics. 



Table 3 Associations between trends 


in emphysema, asthma, and 


pneumonia death rates and dynamics 


of air pollutants in North 


Carolina, 1992-2010 








Potential health-Impact factor 


Emphysema 


Asthma 


Pneumonia 


Ozone, ppb 


0.0061 ±0.0030, P<0.05 


0.0082±0.0056* 


-0.001 1 ±0.00 19* 


Smoking 


0.0493±0.0056t, P<0.000l 


0.0649±0.0l05t, P<0.000l 


0.04l3±0.0034t, P<0.000l 


SOj, ppb 


0.0547±O.OI06t, P<0.000l 


0.0598±0.0l73t, P<0.00l 


0.0309±0.0093t, P<0.00l 


Smoking 


0.0399±0.0074t, P<0.000l 


0.0563±0.0l2|t, P<0.000l 


0.0360±0.0063t, P<0.000l 


NO^, ppb 


0.015310.0062, P<O.OI 


0.0270±0.0094, P<0.005 


0.003010.0053* 


Smoking 


0.0456±0.0090t, P<O.OOOI 


0.051 I±0.0l40t, P<0.00l 


0.0455±0.0076t, P<0.000l 


CO, ppb 


0.0004±0.000|t, P<O.OOOI 


0.0006±0.000|t, P<0.000l 


0.0001 ±0.0001* 


Smoking 


0.0300±0.0083t, P<O.OOI 


0.0349±0.0I29, P<0.0l 


0.0388±0.0074t, P<0.000l 


PM,3, Hg/m^ 


0.015510.0066, P<0.05 


0.01 I6±0.0083* 


0.004410.0063* 


Smoking 


0.04l4±0.0072t, P<0.000l 


0.0329±0.0093t, P<0.00l 


0.0462±0.0067t, P<0.000l 


PM,„, Hg/m' 


0.004510.0039* 


0.0204±0.0058t, P<0.00l 


-0.001510.0035* 


Smoking 


0.0583±0.0069t, P<0.000l 


0.0644±0.0l09t, P<0.000l 


0.0499±0.0057t, P<0.000l 



Notes: For each air pollutant, the effect of smoking was evaluated. The effects of month-to-month fluctuations in disease-specific mortality for emphysema, asthma, and 
pneumonia are not shown in the table, but they also were evaluated for each month. *P>0.05; ^^significant under Bonferroni correction for multiple comparisons. 
Abbreviation: PM, particulate matter. 



618 



submit your manuscript | 
Dovepress 



International Journal of COPD 2014:9 



Dovenress 



Respiratory mortality and changes in air quality 



S? Is 



JS — 
t» 1. 

« IT 



-I .E « 



■o 
c 

3 



3 
M 



O 



o 
Z 



M 
1. 

C 

(4 



>- 



< 
Z 



o 
z 



o" 








z 




Z 












indarc 




£ 

CL 








o 












Z 

JZ 

■t-» 




;cam 




> 


Q. 


T3 

L. 








rt 
■p 






£ 


Star 




effect 


ro 


o" 

Z 


gpm 


o 


O 




p 


c 

s> 


uck; 




s 








:ed 


ndard 


t-dut) 


1977- 


redu< 


eas 


igh 


c 




new 


and 1 


l/J 
P3 




Thei 


cars ; 


For c 


1981 



o 
z 



E 

00 ro 



E 

° o 
- i 

E 
p 



O 

z 



■3 u -p 



0) 00 

E 



O O O O CU 

— Ll_ — Ll_ Ll_ Ll_ -p 



-s 



1^ 



CJ < 



c .y 
o ^ 
U 5 



£ .-^ 



U 



0) 



z ± < 

< — > z 



o 
z 



o 
z 

<- 



^ fj^ Z. 

< > ^ 



-Q . -=- 



"Si) 
o 



Q. _5) ci 



o 



o 

Z 



o 



» 2 C '« c/1 

M 2 o y *l> 

.ii S 4-» 



0) O 



o S 
Z 

E re 
o 



O 



E ^ 



-5 5 O 



O 
z 



O 5! 
E Si 



01 — 0) S 

-Si S - 
"S I - 2 



M O 

t z 



- o 



o 



o s 

Z o 



D "5. H 



1/1 U 

Z Co 
-J3 j2 



re o 

1 i 

O J! 



T3 



CO 



o 



U 

t 

o 
Z 



C — -2! 



o 



OJ 



o 
z 



E 



o 



U o re 



U 



International Journal of COPD 2014:9 



submit your manuscript | 



619 



Dovepress 



Kravchenko et al 



Doveoress 



0) .E 'J 



U a - 



^ M 

o il 



o » 5 
Z a .2 



01 tt 

u. U 



O O 



2 S o 



o 
z 



1 £ 



=5 OS P 



t! 
< 



E 5 



0 c 

(N to 

1 ^ 



t! Si 



0) u 



0) 0) 

> > 

0) 0) 



z" o o~ o 



o" 



o 
u 



o 
z 



o 
z 



E 



O LO 



o 
z 

I i- 



E o 



S! 5 g ^ 
S! ^ £ -S 
&o ^ P F 



5 ■= 



< 
< 



f3 
F 



o" 



0 0-7 

t/) 



E s 
E 

O 



6<j ^ 
E 



ci 2 



u z b; 



U Q 



Z I 



lis 

1- 2 § 

I g s 



s- 5« 0 



< 



5 5^ 0 
-§ I I 

q: Z < 



E Z O 



o 
z 



o o 
o o 



O 

I/) 

6' 



O £ 



< _- 



or 








i3 


E 


S. 


E 


o 














CD 
Q. 




rt 
(U 






X 


ton; 


ires 


per 


o 
o 
o 


requ 


tons 




Iso 


o 
o 


g 




o 




*-> 


o' 




u 


o 


Enerj 


; the ; 


and 1 




i3 


n 


oi 


c 




CM} 




S- 


0 


■u 

(1) 




1^ 




CL 


TD 


«= 


c 


c 

rt 


oa 


g 


year 


u 


o 


.£ 


o 
o 




2 


o' 


pe 


0) 

u 


LO 



§: S E 



8 



U 

c 

o 

z 



< 
z 



I?, 



& s 



.2 £ 



> .2 _L 



< 
> 



< 

3 



Z < 



3 -S 



> 



l_ w 



QJ Pi =f 



rt ™ u 



-c c QJ 



^ -» 



620 submit your manuscript | 

Dovepress 



International Journal of COPD 2014:9 



Dovenress 



Respiratory mortality and changes in air quality 



differences in concentration process (eg, caused by the 
thinning of the air mixing layer in winter), and season-specific 
higher formation of certain compounds, eg, higher nitrate 
formation in the cold season leads to higher levels of NO^ in 
the air.^^ Higher PM levels observed in North Carolina dur- 
ing the summer are of additional concern for health effects 
being exacerbated by hot humid weather, especially during 
heat waves. For respiratory mortality, no threshold effect 
has been identified;^"*'^^ therefore, detailed economic analysis 
is required to evaluate the expenses and benefits of keeping 
the levels of air pollutants extra low. For current regulations 
in the US, it has been shown that control of PM^ ^ emissions 
could result in $100 billion of benefits annually.^'' 

Air quality and emphysema 

In our study, the association between reduced levels of ozone, 
SOj, NOj, CO, and PM^^ and decreased mortality from 
emphysema were observed, with associations for SO^ and CO 
remaining significant under Bonferroni correction. In other 
studies, emphysema outcomes were usually analyzed as a part 
of COPD; nonetheless, our findings on emphysema are in gen- 
eral agreement with these publications. For example, higher 
prevalence of visits to emergency departments for COPD and 
emphysema have been observed for higher SO^ levels" (espe- 
cially among older adults^^); however, some studies showed 
that these associations may be attributable to SO^ serving as 
a surrogate of other substances.^' Few studies are available 
on the effects of outdoor CO on COPD.'"' '" Our results on 
associations between lower CO levels and lower emphysema 
mortality are in agreement with studies that showed increased 
morbidity and mortality risks among patients with COPD.''^^^ 
Note that the impacts of CO could be effectively minimized 
by controlling transportation activities, which accounts for 
more than three-quarters of CO emissions in the US.*^'" While 
in our study associations with PM^ ^ became nonsignificant 
under Bonferroni correction, in other studies higher levels of 
PMj 5 have been associated with higher admissions for COPD 
exacerbation*' and with increased COPD mortality. These 
differences could be due to the fact that the aforementioned 
studies were performed outside the US, had different patterns 
of seasonal fluctuations of PM levels in the air, and also were 
focused on specific populations (ie, older adults). 

Air quality and asthma 

We observed decreasing asthma mortality associated with 
lower levels of NO^, SO^, CO, and PMj^, with the latter 
three pollutants remaining significant under Bonferroni 
correction. These results are in agreement with other studies. 



For example, correlations have been reported between asthma 
mortality and SO^" and NO^^^^^^'' levels, and between asthma 
severity (in children) and CO levels.'"*''" Other studies 
reported that asthma mortality decreased earlier in response 
to improvement of air quality (eg, when compared to emphy- 
sema or chronic bronchitis),^' with a decrease of asthma 
deaths occurring approximately 5 years earlier. 

The effects of PMs on respiratory health and, in particular, 
on asthma have been studied predominantly for associations 
with prevalence of respiratory symptoms^^"*" and emergency 
department visits or hospital admissions.^* '"'"''' It has been 
shown that asthma symptoms were exacerbated even at PMs 
concentrations being 60% below the safety limits for PMs (ie, 
that supposed not to affect the healthy population).'^ However, 
information on associations of asthma mortality with long- 
term exposure to PMs is sparse. In our study, reduction of PMj^, 
(and its seasonal fluctuations) was associated with decreased 
asthma mortality in North Carolina. Previous studies on PM^^j 
showed that elevated levels of PM^^ were correlated with 
hospital admissions for asthma among patients aged 65+ 
years'"^ and children,""'^' and also with increased use of asthma 
medications among patients aged from 8 to 72 years old.'''' 

Because air-quality and asthma-aggravation associations 
are reported from the studies typically performed in a single 
geographic region over a single season, individual study 
results may not be applicable to different populations and to 
longer weather/season cycles.''^ Also, different components of 
PMs (eg, sulfates, nitrates, organic chemicals, metals, and soil 
or dust particles)'^ may have different effects on the respiratory 
system."'-'''''' This makes comparisons between the studies 
challenging and may explain the diversity of results on health 
effects of PMs on both geographic and temporal scales.''* 

Air quality and pneumonia 

In our study, a decrease in pneumonia deaths was associated 
with decreasing SOj levels. Also, when pneumonia was con- 
sidered as the underlying cause of death, lower pneumonia 
death rates were observed for lower CO levels. Some studies 
have linked an acute respiratory disease with higher levels of 
SOj pollution, independently of cigarette smoking,''' while 
later studies have not confirmed these associations (however, 
some results were sensitive to the methods used to estimate 
air-pollutant levels).™" For CO, an association has been 
reported between its increased concentrations and higher 
pneumonia hospitalization.''^ 

While some epidemiological and experimental stud- 
ies have suggested relationships between NO^, ozone, and 
PMs and increased risk for viral respiratory infections,'^ we 
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did not find these associations in our study. Our results are 
in agreement with another study that did not find positive 
associations between PMs and pneumonia deaths (they found 
associations only for the group of never-smokers).'' However, 
most of the studies were performed on pneumonia morbidity 
(including hospitalizations and emergency department visits), 
while our study was on mortality. Also, multiple reports on 
associations between pneumonia risk and PMs levels come 
from international studies, eg, from Europe (where PMs 
levels peak in winter), while on the East Coast of the US they 
typically peak in summer,'" as we also observed in our study. 
While pneumonia is more frequent in late fall and winter, 
the relationships between outdoor air quality and health are 
supposed to be stronger in summer, when people spend more 
time outdoors. A study from Boston also supports our find- 
ings; no associations with pneumonia hospital admissions 
were found in summer, while in winter the largest effect on 
pneumonia morbidity was reported not for PMs but for black 
carbon (a surrogate for traffic particles: 14.3% increase of 
pneumonia hospitalizations for 1.7 |J,g/m^ increase of black 
carbon).''^ Higher risk of morbidity and mortality from acute 
respiratory infections has been also reported for children 
exposed to PM^jj.^^'^"**' In our study, we did not estimate 
mortality risks specifically for children; future studies will 
be performed for age-groups that are potentially at highest 
risk (ie, children and older adults). 

Methodological aspects 
and study limitations 

In our approach, the number of observations sufficient to 
estimate model parameters was achieved by incorporating 
monthly changes of air-pollutant levels and respiratory 
mortality. One advantage of this approach is that the 
unobserved heterogeneity due to other factors (such as 
socioeconomic status, quality of health care, migration) is 
minimal, because these factors do not essentially vary from 
month to month. In contrast, this unobserved heterogeneity 
is typical for ecological studies with area-based design, and 
could result in the occurrence of additional biases if these 
variables are not sufficiently controlled. 

One example of such a factor is the time trend describ- 
ing improvements in the treatment of respiratory disease 
that occurred during the recent two decades and which 
contributed to decreasing trends of mortality from emphy- 
sema, asthma, and pneumonia. Both improved air quality 
and vaccinations against pneumonia could lead to fewer 
hospital admissions, eg, pneumonia age-adjusted death 
rates started declining in the late 1990s, while the hospital 



discharge rate did not change significantly for patients older 
than 1 5 years. Although our approach with measurements 
at the month level minimizes the bias from this time trend 
(because only a 12th of our measurements reflect the annual 
time trend), improvements in treatment (as well as factors 
other than air pollution and smoking with significant time 
trends) should be taken into account in further studies. 
For example, further analysis of disease-specific visits to 
emergency departments would be important to validate 
the role of improved medical care in observed respiratory 
disease trends. 

Other factors, such as changes in socioeconomic status, 
can also impact the dynamics of disease-specific mortality 
rates. However, it has been reported that for social factors, 
as well as for race, the effects of modification, eg, of PMs 
(ie, PMjj,) on total mortality were weak.^" 

In our study, the time pattern of smoking was chosen to 
reflect annual trends in respiratory mortality in addition to air 
pollution. Inclusion of one additional variable measured annu- 
ally (ie, not on a monthly basis) could result in difficulty in 
distinguishing the effect of this variable and smoking. Smoking 
was chosen because its patterns are concordant with patterns of 
respiratory mortality, and because of many substantive results 
on the role of smoking in respiratory mortality (eg, findings that 
both smoking and exposure to air pollutants [eg, PM^ J could 
exacerbate respiratory diseases). ^^''^ In our study, smoking 
had a significant stable effect on the dynamics of respiratory 
mortality from all three studied diseases. However, it can also 
reflect possible impacts of other variables with similar to smok- 
ing time trends and associations with respiratory mortality. 
Better evaluation of smoking effects (including synergistic 
effects of smoking and air pollutants) could be achieved in 
studies with individual records on smoking status. 

Study designs based on individual measurements of envi- 
ronmental exposure and health outcomes (which are classic 
epidemiologic approaches) would be helpful for improve- 
ment of the quality of estimates. However, such approaches 
are expensive and complex, in part due to the difficulty of 
measuring subjects' exposure to the relatively low levels 
of pollutants in the air. Some studies on the use of outdoor 
monitoring-station data (compared with the personal indoor/ 
outdoor-exposure monitors) demonstrated that personal 
exposure to pollutants of outdoor origin was more closely 
related to outdoor air-pollutant levels than interpretations of 
personal monitoring data.^*-^' Furthermore, the frequently 
high correlations between levels of certain pollutants in the 
air also make it difficult to identify the impact of a single 
agent on human health." 
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Changes in diagnostic criteria of respiratory diseases 
that happened during last two decades primarily affected 
the trends of disease incidence; however, in part, mortal- 
ity trends were also affected. In children, diagnoses can 
transfer from chronic bronchitis and pneumonia to asthma, 
thus contributing to increasing trends in asthma prevalence 
(with its recent stabilization) and health care utilization.'" If 
the person dies from pneumonia, but also had an underlying 
condition of which the pneumonia was probably a result, than 
that underlying disease but not pneumonia is considered the 
cause of death in the death certificate, and thus fewer deaths 
are directly attributable to pneumonia.*"" Although asthma 
death rates increased from 1980 to the mid-1990s, replaced 
ICD codes from the ninth to the tenth revision makes it chal- 
lenging to evaluate the decline in asthma mortality since the 
late- 1 990s." With regard to this problem, it has been shown 
that decline in asthma mortality that occurred from 1998 to 
1999 included approximately 11% of decline that resulted 
from the changes during the ICD codes transition; then, 
under ICD- 10, asthma death rates continued declining." 
Because no definitive asthma laboratory tests exist, asthma 
estimates rely on the physician, who also should accurately 
attribute the cause of death to asthma; therefore, the reli- 
ability of the death certificates has been questioned (eg, 
for the chance of misreporting the cause of death in older 
persons with comorbid conditions). Large well-designed 
studies have concluded that asthma death coding has 99% 
specificity and low sensitivity (42%), and asthma as a cause 
of death was underreported in preference to COPD in all 
age-groups."-'^ 

Conclusion 

We observed temporal regional associations between 
long-term dynamics of decreasing death rates of emphy- 
sema, asthma, and pneumonia and reductions of the levels 
of certain air pollutants in North Carolina. Our results 
support the hypothesis that improvement in air quality, 
especially declines in SO^, CO, and PM^u levels in the air, 
contributed to the improved respiratory health of the North 
Carolina population. Since other factors (in addition to the 
studied air pollutants) might also account for improved 
health outcomes, ultimately caution should be exercised 
in inferring cause-effect relations. 
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Table SI Results of the sensitivity analysis 



Potential health-Impact factor 


Emphysema 


Asthma 


Pneumonia 




Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


Ozone, ppb 


0.005610.0029, P>O.OS 


0.005210.0050, P>0.05 


-0.001 110.0019, P>0.05 


Smoking 


0.043310.0057, P<O.OOOI 


0.042710.0097, P<0.000l 


0.041 110.0036, P<0.000l 




Analysis 2 (summer) 


Analysis 2 (summer) 


Analysis 2 (summer) 


Ozone, ppb 


0.004310.0069, P>0.05 


-0.000410.0085, P>0.05 


-0.004610.0024, P>0.05 


Smoking 


0.077710.0159, P<O.OOOI 


0.1 14010.0205, P<0.000l 


0.040310.0051, P<O.OOOI 




Analysis 3 (winter) 


Analysis 3 (winter) 


Analysis 3 (winter) 


Ozone, ppb 


0.009210.0037, P<O.OI 


0.000410.0098, P>0.05 


0.005210.0056, P>0.05 


Smoking 


0.037710.0066, P<O.OOOI 


0.010410.0153, P>0.05 


0.024 1 10.0 1 0 1, P<0.05 




Analysis 4 (underlying) 


Analysis 4 (underlying) 


Analysis 4 (underlying) 


Ozone, ppb 


0.003910.0031, P>0.05 


-0.000310.0084, P>0.05 


0.005410.0049, P>0.05 


Smoking 


0.052410.0056, P<O.OOOI 


0.068210.0156, P<0.000l 


0.081410.0097, P<0.000l 




Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


SOj, ppb 


0.050210.0108, P<O.OOOI 


0.028910.0159, P>0.05 


0.033110.0094, P<0.00l 


Smoking 


0.036110.0075, P<O.OOOI 


0.037510.0109, P<O.OOI 


0.037910.0064, P<O.OOOI 




Analysis 2 (summer) 


Analysis 2 (summer) 


Analysis 2 (summer) 


SOj, ppb 


0.055110.0189, P<0.05 


0.053510.0231, P<0.05 


0.002710.0077, P>0.05 


Smoking 


0.065110.0145, P<O.OOOI 


0.097810.0188, P<0.000l 


0.034510.0052, P<0.000l 




Analysis 3 (winter) 


Analysis 3 (winter) 


Analysis 3 (winter) 


SOj, ppb 


0.082310.0243, P<O.OOI 


0.029810.0357, P>0.05 


0.059610.0250, P<0.05 


Smoking 


0.038610.0167, P<O.OOOI 


0.043210.0247, P>0.05 


0.019410.0170, P>0.05 




Analysis 4 (underlying) 


Analysis 4 (underlying) 


Analysis 4 (underlying) 


SOj, ppb 


0.035810.0126, P<0.005 


0.038710.0259, P>0.05 


0.109410.0193, P<0.000l 


Smoking 


0.047110.0087, P<O.OOOI 


0.072410.0180, P<O.OOOI 


0.066610.0140, P<0.000l 




Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


NOj, ppb 


0.015910.0062, P<O.OI 


0.028110.0084, P<O.OOI 


0.002910.0053, P>0.05 


Smoking 


0.036710.0094, P<O.OOOI 


0.017910.0129, P>0.05 


0.046910.0079, P<0.000l 




Analysis 2 (summer) 


Analysis 2 (summer) 


Analysis 2 (summer) 


NOj, ppb 


0.043410.01 14, P<O.OOOI 


0.016010.0163, P>0.05 


0.004910.0051, P>0.05 


Smoking 


0.041010.0167, P<O.OI 


0.098710.0235, P<0.000l 


0.030610.0067, P<0.000l 




Analysis 3 (winter) 


Analysis 3 (winter) 


Analysis 3 (winter) 


NO^, ppb 


0.013510.0149, P>0.05 


0.022410.0196, P>0.05 


-0.01 10110.0145, P>0.05 


Smoking 


0.062110.0203, P<0.005 


0.035210.0263, P>0.05 


0.059110.0199, P<0.005 




Analysis 4 (underlying) 


Analysis 4 (underlying) 


Analysis 4 (underlying) 


NO^, ppb 


0.002410.0072, P>0.05 


0.038510.0133, P<0.005 


0.009810.01 14, P>0.05 


Smoking 


0.059010.0104, P<O.OOOI 


0.047710.0204, P<0.05 


0.098610.0171, P<0.000l 




Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


CO, ppb 


0.000410.0001, P<O.OOOI 


0.000210.0001, P>0.05 


0.000210.0001, P<0.05 


Smoking 


0.029910.0084, P<O.OOI 


0.034210.0120, P<O.OI 


0.038610.0073, P<0.000l 




Analysis 2 (summer) 


Analysis 2 (summer) 


Analysis 2 (summer) 


CO, ppb 


0.001310.0003, P<O.OOOI 


0.001710.0004, P<O.OOOI 


0.000210.0001, P>0.05 


Smoking 


0.051010.0149, P<O.OOI 


0.076610.0179, P<O.OOOI 


0.029810.0057, P<O.OOOI 




Analysis 3 (winter) 


Analysis 3 (winter) 


Analysis 3 (winter) 


CO, ppb 


0.000510.0001, P<O.OOI 


0.000710.0002, P<O.OOOI 


0.000110.0002, P>0.05 


Smoking 


0.026710.0202, P>0.05 


-0.017310.0258, P>0.05 


0.045410.021 1, P<0.05 




Analysis 4 (underlying) 


Analysis 4 (underlying) 


Analysis 4 (underlying) 


CO, ppb 


0.000 1 10.000 1, P>0.05 


0.000810.0002, P<0.000l 


0.00 1 010.000 1, P<0.000l 


Smoking 


0.050 1 10.0 1 0 1, P<O.OOOI 


0.033410.0194, P>0.05 


0.035210.0150, P<0.05 




Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


PM,3, ng/m^ 


0.015510.0066, P<0.05 


0.01 1610.0083, P>0.05 


0.004410.0063, P>0.05 


Smoking 


0.041410.0072, P<O.OOOI 


0.032910.0093, P<O.OOI 


0.046210.0067, P<0.000l 




Analysis 2 (summer) 


Analysis 2 (summer) 


Analysis 2 (summer) 


PM^ j, ng/m^ 


0.020710.01 13, P>0.05 


0.001410.0105, P<0.05 


0.001610.0045, P>0.05 


Smoking 


0.057810.0152, P<O.OOOI 


0.079710.0141, P<O.OOOI 


0.035210.0060, P<0.000l 


(Continued) 
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Table S I (Continued) 



Potential health-impact factor 


Emphysema 


Asthma 


Pneumonia 




Analysis 3 (winter) 


Analysis 3 (winter) 


Analysis 3 (winter) 


PM,,, ng/m' 


0.022410.0168, P>0.05 


0.019510.0186, P>0.05 


0.003910.0214, P>0.05 


Smoking 


0.049810.0154, P<0.00l 


0.004910.0168, P>0.05 


0.044610.0193, P<0.05 




Analysis 4 (underlying) 


Analysis 4 (underlying) 


Analysis 4 (underlying) 


PMjj, ng/m' 


0.003010.0072, P>O.OS 


0.003010.0137, P>0.05 


-0.004710.0101, P>0.05 


Smoking 


0.056010.0078, P<O.OOOI 


0.036610.0154, P<0.05 


0.066710.0106, P<0.000l 




Analysis 1 (ICD-9/I0) 


Analysis 1 (ICD-9/10) 


Analysis 1 (ICD-9/10) 


PM,„, Hg/m= 


0.002510.0039, P>0.05 


0.012510.0053, P<0.05 


-0.001210.0035, P>0.05 


Smoking 


0.052110.0072, P<0.000l 


0.039510.0101, P<0.000l 


0.050810.0059, P<0.000l 




Analysis 2 (summer) 


Analysis 2 (summer) 


Analysis 2 (summer) 


PM,„, ng/m' 


0.016910.0071, P<0.05 


0.024310.0083, P<0.05 


-0.002910.0026, P>0.05 


Smoking 


0.071410.0142, P<0.000l 


0.102010.0176, P<0.000l 


0.037710.0050, P<0.000l 




Analysis 3 (winter) 


Analysis 3 (winter) 


Analysis 3 (winter) 


PM,„, Hg/m' 


-0.010410.0143, P>0.05 


0.040710.0180, P<0.05 


-0.013510.0143, P>0.05 


Smoking 


0.082810.0170, P<0.000l 


0.029710.0213, P>0.05 


0.056910.0162, P<0.0005 




Analysis 4 (underlying) 


Analysis 4 (underlying) 


Analysis 4 (underlying) 


PM,„, Hg/m' 


-0.006310.0047, P>0.05 


0.025610.0084, P<0.005 


0.014010.0074, P>0.05 


Smoking 


0.066210.0078, P<0.000l 


0.067910.0163, P<0.000l 


0.099310.0133, P<0.000l 



Notes: The following factors were tested; the potential effect of International Classiftcatjon of Diseases (ICD) code changes (from ICD-9 to ICD- 1 0) (analysis I ), the effects of 
air pollutants on mortality during the summer (analysis 2) and winter (analysis 3), and the association when only underlying causes of death contributed to the cause-specific 
death rates (analysis 4). 
Abbreviation: PM, particulate matter. 
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